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We discuss a general approach to a realistic theory of the electronic structure in materials 
containing correlated d- or f- electrons. The main feature of this approach is the taking into account 
the energy dependence of the electron self-energy with the momentum dependence being neglected 
(local approximation) . It allows us to consider such correlation effects as the "non- Fermi-step" form 
of the distribution function, the enhancement of the effective mass including "Kondo resonances", 
the appearance of the satellites in the electron spectra etc. To specify the form of the self-energy, it 
is useful to distinguish (according to the ratio of the on-site Coulomb energy U to the bandwidth W) 
three regimes - strong, moderate and weak correlations. In the case of strong interactions (U/W > 1 
- "rare-earth system") the Hubbard-I approach is the most suitable. Starting from an exact atomic 
Green function with the constrained density matrix n„„/ the band structure problem is formulated 
as the functional problem on n^^/ for f-electrons and the standard LDA-functional for delocalized 
electrons. In the case of moderate correlations (U/W~l -"metal -insulator regime, Kondo systems") 
we start from the d= oo dynamical mean field iterative perturbation scheme (IPS) of G. Kotliar et. 
al. and also make use of our multiband atomic Green function for constrained Um^/. Finally for 
the weak interactions (U/W < 1 -"transition metals") the self-consistent diagrammatic fluctuation- 
exchange (FLEX)-approach of N. Bickers and D. Scalapino is generalized to the realistic multiband 
case. We presents two-band, two-dimensional model calculations for all three regimes. A realistic 
calculation in Hubbard-I scheme with the exact solution of the on-site multielectron problem for 
f(d)- shells was performed for mixed-valence 4f compound TmSe, and for the classical Mott insulator 
NiO. 
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I. INTRODUCTION 



A general accurate description of the electronic structure of materials with correlated electrons has yet to be 
developed. Such materials include the high-Tc and colossal magnetoresistance (CMR) materials, as well as the 
mixed-valence and heavy-fermion compounds. All these systems demonstrate essentially many-particle (correlation) 
features in their excitation spectrum and ground-state properties, the usual language of one-electron band theory 
being inadequate to describe such features even qualitatively: e.g., the problem of Mott insulators, the heavy-fermion 
behavior in some rare-earth compounds, satellites and "mid-gap states" in electron spectra etc. (see e.g. recent reviews 
[Q-Q). Such effects as the metal-insulator transition, Kondo effect, and others, which helps to understand the basic 
physics in these strongly correlated materials, is usually considered in the framework of simplified models such as the 
Hubbard model, Anderson model, s-f exchange model and other correlation models. Nevertheless the complexity of the 
crystals containing 10-15 different atoms per unit cell, interactions between electronic and lattice degrees of freedom 
demands a more detailed investigation of the energy bands in such systems. The only general first-principles approach 
that takes into account in practice specific peculiarities of the electronic structure in real compounds are those based 
on the density functional theory (DFT) [^|. The vast majority of practical DFT applications today are based on the 
local density approximation (LDA), which treats the exchange-correlation (XC) part of an effective single-particle 
DFT potential as a density dependent XC-potential, taken from the exact quantum Monte-Carlo (QMC) results for 
the homogeneous electron gas. There are many successes but also some failures of LDA approach related to 
the simple fact that in cases where some portion of electronic structure is better described in terms of atomic-like 
electronic states, the homogeneous electron gas approximation is not a good starting point. Another limitation of 
LDA theory is that it is only a ground state scheme and the one-particle band-structure itself has, generally speaking, 
no proper meaning. Recently the time-dependent (TD-LDA) approach has been applied for calculations of excitation 
energies but the TD-LDA effective potential is not known as good as LDA one. On the other hand, exact 

QMC calculations for real materials which have tedious first-principle Hartree-Fock band structure as zero order 
approximations is still a challenging problem of solid state theory M- . 



In this situation it is useful to have simple and accurate scheme that could still capture the most important properties 
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of real electronic structure and at the same time could take into account the most important correlation effects. One 
of the first successful approach in this line was the GW-approximation for quasiparticle spectra in solids with 
the self-energy related to a "bare" Green function (G) and a screened Coulomb interaction (W). Self-consistent GW 
calculations basing on the LDA band structure gives much better description of a Mott insulator such as NiO than does 
pure LDA Still, the non-local Coulomb interactions make such type of calculations really time consuming. For 
the purpose of only band-structure investigations one could use a simplified time-independent GW-scheme or so-called 
screen-exchange (sX-LDA) approach . The later approach has of the same drawbacks as does Hartree-Fock 

approximation and did not suitable for strongly correlated systems. A different way to incorporate some correlation 
effect in the systems with localized d- or f-states was successfully done in so-called LDA-I-U method ||l6| . In this case 
a simple mean-field Hubbard-like term is added to LDA functionals for the localized state and care must be taken for 
correction of the LDA double-counting This approach can also be viewed as a density functional theory, since 
the U terms that depend on occupation number for localized electrons is a functional of the total density. So one just 
uses the LDA- functional for delocalized electrons and improved LDA-I-U functional for localized atomic- like states. 
This approach produces a more reliable description of the electronic and crystal structure of correlated materials with 
charge, spin, and orbital ordering than does the LSDA scheme But LDA-I-U scheme, as well as the approach 

based on the so called self-interaction corrections (SIC) |l^, has one intrinsic shortcoming related with a mean- 
field approximation. It is well known (see e.g. [^0(12 1| ]) that the most interesting correlation effects in quasiparticle 
spectra such as the mass enhancement, damping, the difference of the distribution function from the "Fermi step" are 
connected with the energy dependence of the self-energy S(w), so one needs to generalize LDA+U approach to include 
dynamical effects. Such scheme we would like to call "LDA++" One can mentions few successful attempts in 

this directions: quasiparticle (QP) band structure calculation of Fe, Co and Ni p3|-p5t as well as heavy fermion system 
p6| , p7[ using simplest second order local approximation for the self-energy; QP-band structure of NiO ||2^,^ using 
3-body Faddeev approximation; RPA-like approach for HTSC |Q and non-crossing approximations (NCAjfor Kondo 
systems pl| ]. At the same time, a criterion for the applicability of specific approximations used in these works was 
not clear. 

In this paper we propose a general scheme for LDA-f-|- band structure calculations for real materials with the 
different strength of electronic correlations. It is not very efficiently from the computational point of view (as well as 
not very reasonable from the purely theoretical one) to use the only LDA++ scheme for materials with a different 
electron-electron interactions. In accordance with the ratio of average on-site Coulomb parameter U to the relevant 
valence band width W, it is useful to distinguish three regimes of weak, moderate and strong correlations. 

The simplest case of weak correlations (U/W<1 - "transition metals") we could use the self-consistent diagrammatic 
approach. The most convenient way is the conserving fluctuation exchange (FLEX) approximation and we will use 
the multiband generalizations for LDAh — h weak correlation scheme. The characteristic feature of this renormalized 
band regime that no additional states appears in the electronic structure due to interactions or, more exactly, there 
is one-to-one correspondence of quasiparticle states with and without interaction. Roughly speaking, the shape of 
energy bands may be changed but there are no band splitting, or presence of additional bands. The most interesting 
physical phenomena in this case are the renormalization of the effective masses, flattening of Van Hove singularities 
etc. 

In the case of very strong interactions (U/W>1 - "rare-earth system") we will start with exact atomic Greens 
function for f-states and use the Hubbard-I approximation (HIA) to analyze the spectrum of f-systems. This 
approach also may be applied to such d-systems as Mott insulators with very narrow d-bands. For this situation the 
electronic structure of solids will combine the many-body structure of f(d)-ions and broad bands from delocalized 
electrons. In this case such phenomena as multiplet structure, satellites in photoelectron spectra, the narrowing of 
the electron bands depending on the magnetic ordering etc are the subjects of main interests. 

In the most difhcult case of strongly correlated physics (U/W ~1 - "metal- insulator transition" regime or "Kondo 
systems") we will use the interpolation scheme based on dynamical mean held theory (DMFT) In this situation 
we will have a "three-peak structure" from a single correlated band, consisting of upper and lower Hubbard bands 
and Kondo resonance near the Fermi energy. Such scheme is the most accurate, but also the most time consuming, 
and it is difficult to make a self-consistent calculations for a large system. The other point of view of the different 
LDA++ schemes could be related to the different energy scales for the spectrum of correlated materials: if one is 
interested in the large energy scale that HIA-approximation is sufficient for spectroscopic purposes. If we like to 
describe low-energy scale of a system like doped Mott insulators or mixed valence systems then DMFT approach is 
the most adequate. 

The common feature of all LDAH — h methods is the matrix form of self-energy since electron-electron correlations 
can not be diagonalized neither in band index n nor in orbital indices Im. This peculiarity of multi-band Hubbard 
interactions are normally ignored and only few examples of matrix self-energy exist for transition metal with LDA 
second-order perturbation scheme |23|j2^ ], transition-metal oxides within 3-body Faddeev approximation [ p9| and for 
two-band Hubbard model for investigation of orbital and magnetic instabilities Q . 
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This paper is organized as follows. In Sec. II. we will give a general description of different correlations scheme to 
the dynamical mean field band structure calculations. The simple two-bands model (Sec. Ill) will illustrate in practice 
all LDA-I--I- methods. In the Sec. IV we will give an example of first-principle calculations of mixed-valence system 
TmSe and classical Mott insulator NiO within Hubbard- 1 approximation to LDAh — h. Finally we summarize our 
results in the sec V. 



II. LDA-I-+ METHODS 



The "Kohn-Sham" energies of one-particle LDA states cannot be considered as the quasiparticle energies in the 
sense of many-particle theory (see e.g. [ p0|j2l| ). In the LDA-I--I- approach they considered only as the "bare" energies 
and are supposed to be renormalized by the correlation effects. Of course they contain already some part of the 
correlation effects but only those which may be considered in the local density approximation. The most important 
"rest" in strongly correlated system is the correlations of the Hubbard type |Q due to the intra-site Coulomb repulsion. 
Therefore our starting point is the same as in the LDA+U approach. We proceed with the Hamiltonian 

H = f^^ r+ r- TP , , r+ r+ r , ,r > (^] 

ija{m} zjcrm} 

where the represents different crystal sites, {m} label different orbitals and the {a} are spin indices. Coulomb 
matrix elements are defined in the usual way: 



J ^'"^^ ('■'^^-(^ - (2) 

here Vee(r — r') is the screened Coulomb interactions and iprn,{j^) are localized on-site basis functions (the site index 
being suppressed). 

In this case all the orbitals are assumed to belong to the "correlated" set, while in real materials like high-Tc 
compounds (e.g., YBa2Cu307 ) we may define as a first approximation only the 3d orbitals as correlated ones. 
Therefore it is more reasonable to rewrite Eq.|l| in the form of LDA-I-U Hamiltonian: 

vr^ljLDA,}_ \^ TTi r+ r+ r , ,r , (V\ 

dc ^ 2 ^ ^mimam'jmj imicr imaf ""2'^ ^h*^ ^ ' 

i{(7m} 

Note that index i in the second sum of Eqj^ running only for "correlated" sites, and orbital indices {m} only for 
"correlated" states (e.g. 3d or 4f) while the first LDA-term: 

ttLDA _ ,Z3 + _ p 

^dc — / , "'mim2a'^imia'^]rn2a J^dc 

ijcr{m} 

contains all sites and orbitals in the unit cell. Here K^_^^^^ are the "one-particle" Hamiltonian parameters in the 
(spin-polarized) LDA, Edc is the double counting correction for average Coulomb interactions in L(S)DA [ pTf : 

Edc = ]^Und{l - rid) - ]^J[nd]{l - rid]) + "-<ii(l - "-di) 

with U and J being average Coulomb and exchange interactions and rid = rid] + ridi is the total number of correlated 
d(f)-electrons. 

One non-trivial problem is to find an efhcient way to compute the one electron Hamiltonian h^^^ in a minimal 
orthogonal basis set. Orthogonality of the basis functions is required for the use of second quantization form of the 
effective Hamiltonian (Eq.|l|). Since many-body {U) part of the problem is an order of magnitude more time consuming 
than LDA one, we need to use a minimal basis set and integrated out all high-energy degrees of freedom (out of the 
±U range). One of the best LDA-methods for such scheme is the Linear Muffin-Tin Orbital (LMTO) method [ p5[ , 
which could give the orthogonal down-folded one-electron tight-binding Hamiltonian. In this case LDA calculations 
corrected for double counting, produce the first-principle hopping tij in many body Hamiltonian. 

Now we can describe methods to efficient calculations of quasi-particle (QP) spectra for the LDA-I-U Hamiltonian. 
In this sense our approach is not any more density functional theory and one could benefit from the possibility to use 
the information on QP-band structure comparing with different "excitation" experiments. 
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A. Multi-band FLEX 



In this section we generalize the FLEX equations ||3^ for the purpose of multiband LDA++ scheme. We will not take 
into account a momentum q-dependence of the self-energy, although in the FLEX approximation it is straightforward 
to include it in all the following formulas. The numerical computation of the (q, u;)-dependent self-energy is time- 
consuming in multiband case [3^j3^ . To unify the approximations for all our LDA++ schemes we will not include 
explicitly the q-dependence in the FLEX formalism. 

First of all one needs to symmetrize the bare vertex matrix U over different fluctuation channels: particle-hole 
(density - U'' and magnetic - U™) and particle-particle (singlet - U'* and triplet - U*) vertex matrix: 

, , ^2U ' ' - U 

mim-^m2m^ mim2m^m^ ■mim2rn^m-^ 

mim^m2m^ 2 'm,\m^m^m2 m\m^rn2m^l 

The one-electron Green-function is defined through the following equation: 

here is chemical potential, a;„ = (2n-f 1)//? are Matsubara frequencies and /3 = Xjk^T is inverse temperature. The 
frequency independent Hartree-Fock part is 

^mm'cr ^ ^ i^mniim' m2 ^ ^ ^m\'m2 ^mmim2m' ^mim2 ) (^) 
mim2 cr' 

and corresponds to the rotationally invariant LDA+U method |Q . 

It is useful to write the multi-band FLEX equations using matrix-vector notation for different Coulomb matrix 
vertices and vector Green function. We will use a combine index: a — {m, m } and defined the vector Green function 
as well as matrix interactions in the following way: 

For simplicity we first write equations for non-polarized spin states and omit the spin indices. In this case the Hartree- 
Fock approximation Eq.^ can be rewritten in the form of a matrix-vector product only with "density" Coulomb 
interaction: 



where the occupation matrix is defined as: 



Using the single-site Hubbard interactions one obtains a "local form" of FLEX equations in the frequency [uo) - time 
(r) space. It is very efficient to use fast- Fourier transforms with periodic boundary condition p^ . Time-frequency 
spaces are connected by 

r/5 



G(iw„) = / e''^"^G(T)dT 
Jo 



We will try to keep this dual {lo — r) notation to stress the numerical implementation of this LDA-|--|- scheme. We 
write the approximation for self-energy in the "G'VF"-like form : 

S(t) = I?(r) * G(t) (5) 
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where symmetrized fluctuation Ty(ici;)-potential defined as: 



and total fluctuation potential consists of the second-order term, as well as particle-hole and particle-particle contri- 
butions: 

V{iu!) = V2{iLu) + Vph{iuj) - Vpp{-iuj) 

All these contributions can be expressed in terms of bare {Dq, Mq, So,To) and renormalized {D,M,S,T) channel 
propagators. The second order potential for the non-magnetic case is 

V2{iuj) ^ U * Doiiuj) * U'^ (6) 

while the particle- hole potential is expressed through the density and magnetic fluctuations: 

Vph{iuj) = ^U'^ * [D{ioj) - Do{iuj)] * U'^ + yj"' * [M{iuj) ~ Mo(iw)] * t/"* 
Finally the particle-particle contribution to fluctuation-exchange potential is: 

%p{iuj) = t/" * - * + 3C/* * [f (iw) - fo(iw)] * [/* 

If one defines the particle-hole (x) and particle-particle (tt) "empty loop" susceptibilities: 



we can write with this notations for susceptibilities the bare channel propagator matrices in the following form for 
the density and magnetic part: 



^0 = Mo = X 



and for singlet and triplet bare propagators: 



The total channel propagators (Ra where A={d, m, s, t}) have to be found from the RPA-like matrix inversion: 

Rx{iuj) = [T + Rl{iuj) *U^]'^ * Rl{i^^) 

The derivation of the complete expression for the FLEX self energy for spin-polarized case with taking into account 
all the channels are rather cumbersome. Since we will not use here these complicated expressions, they will be 
discussed in details elsewhere 

B. Hubbard-I approximation 

Historically Hubbard-I approximation p3| was the first and the simplest approximations for a strongly correlated 
one-band model. It has, however, many inconsistencies (see e.g. discussion in [Q). For example it is not conserving 
(the self energy cannot be represented as a functional derivative of the generating functional with respect to the 
Green function) and therefore does not obey the Luttinger theorem and other "exact" Fermi-liquid properties. For 
the half-flUcd nondegenerate Hubbard model it always gives a gap in the energy spectrum, even for small U. This 
means that HIA is completely inapplicable for small and medium interactions. But at the same time it gives a correct 
picture of the electron spectrum in the narrow-band limit. Therefore it seems to be very useful in 4f systems with a 
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very strong degree of localization of the electron states. Applying to some real systems in the framework of LDA++ 
approach, HIA-scheme could give (as it will be shown below) an effective and non-trivial descriptions of many-body 
multiplet effects. 

To introduce Hubbard-I type approximation in degenerate case it is convenient to exploit the so called atomic 
representation and Hubbard X-operators (see [^3|j|,Q) 

where ii,v are multielectron states of the site i as a whole (configuration and multiplet indices). In terms of X- 
operators the atomic Hamiltonian has very simple form 

On the other hand, the inter-site transfer Hamiltonian, which has very simple (bilinear) structure in terms of the 
operators c^/^, Cm^, also can be expressed in terms of X-operators by the relations c,nci = (a^ \cma\ v) X^^ and 
similarly for c^'^- In the limit of a very strong interaction it is convenient to calculate the Green function via X- 
operators (using the decoupling procedure or special diagram technique for AT-operators pO|) and then transform 
to electron operators. HIA corresponds to the following expression [ p3|j44| ]: 



G-^ {iLj) = [G°*(zc^)] ^ -t 



where t is the matrix of transfer integrals. In the limit of very small t this expression describes the arising of separate 
bands from each intraatomic transition with the change of the electron number from unity. It is the picture that 
seems reasonable for e.g. rare-earth materials with very narrow 4f-band. The bands always appear to be narrowed. 
Indeed, if in the vicinity of the pole iuj = Eq the atomic Green function can be represented as 



J - £0 

the effective transfer Hamiltonian for this "Hubbard band" will be Zgt instead of t. 

In terms of LDA-f -I- multiband approach HIA for Green function has the following form 

G7^ . , (iuj) = \(iuj + (1)6 ' - E"* , (iuj)] 5,, - h'^ , (8) 

im,jm ' '1^1 mm mm /J 'J mm a ^ ' 

To obtain this Green function, we need to solve by an exact diagonalization (ED) technique the atomic many-electron 
problem: 

>= Ef\v > 

with the effective atomic Hamiltonian for d- or f-states: 

fjat ^ + I 1 r r , ,c'^ c'^ ,c > ,c ' (9) 

/ y ^mm ma m cr ' 2 / y mim2m-^m^ im-ia^^jy^^^' ^im^a im^a \ y 

mm' a {am} 

Here is the matrix of atomic energies which in principle can include non-diagonal terms. The latter is naturally 

come from LMTO-TB effective Hamiltonian which has a diagonal part of h^i„ as a result of transformation to an 

rl — I. n mm a 

orthogonal basis set [|35| . Diagonalization of atomic Hamiltonian Eqj9|is not a big problem for a standard workstation, 
since it is equivalent to 5- and 7-site Hubbard model in ED-scheme Q for d- and f-states. 

Using eigenfunctions and eigenvectors of Hamiltonian (EqJ^, the exact atomic Green-function can be found by the 
standard definition [Ool: 



1 \Cma \ v) Lv , LL j 

, (,^) ^ i y ' ^\ ^/ ^e-P^. + e-^^") (10) 

mma^ > iuj + En ~ ^ ' ^ ' 

where Z = E^.e^'^^-. 

Finally which is needed for HIA-approximation is found from the following expression: 
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E"* , {iuj) = iLjS^^> - e< - (G"'*)-i , (iw) (11) 

mm cr ^ ^ vara rara \ ^ mm ^ ^ ' 

Now the HIA-approach to LDAh — h may be formulated as a functional for atomic density matrix: n^^i with a 
constraint (for e^^> ) 

"mm—^E^mrn'M + ^^mm' (12) 

(see ||20|) having the same density matrix for d- or f-electrons as in the crystal as for the corresponding site and 

orbital element of the Green function, Eqj8| 



C. DMFT-multiband scheme. 



A great success of dynamical mean field (d=oo) approach to the theory of correlated systems |^ shows that 
probably this scheme can be the most accurate for the calculations of the self energy from local description of electron 
fluctuations, at least in the vicinity of the metal-insulator transition. We will use this scheme for real crystals as a best 
local approximation. The DMFT scheme is based on the "cavity method" or the solution of the effective impurity 
problem, which corresponds to subtraction of the local self energy only on the one atom in question. In appendix 
A we show the equivalence of the cavity and impurity methods for matrix multiband Hamiltonians. It was realized 
recently that the success of DMFT in the one-band half-filled Hubbard model with simplest second-order self-energy 
is related to the fact that both small and large U-limits are exact in this case This is not true for non-integer 
filling or for the multiband case. The elegant iterative perturbation scheme (IPS) for non-integer one-band Hubbard 
model was proposed recently |4^] and gives almost perfect agreement with ED- and QMC-results. For the case of 
multi-band with non-integer occupations the problem is much more severe and the existing IPS-generalization |4^ ] 
does not produce good results for large doping. Here we use the main idea of the original IPS-method ]4^ ] and propose 
another version of multiband DMFT which is based more on the numerical solutions of corresponding atomic problem 
than the approximate analytical one used in . 

The impurity problem for the "bath" Green function reads: 



where the local Green function is defined through Brillouin zone sum 

Gmm- i}^) ^ ^^G. 



(k, ILO) 



here Nk is the total number of k-points. Alternatively one may perform the k-integration using a complex-tetrahedron 
scheme . We introduce here according to the " local" impurity chemical potential /iq to satisfy the condition 



+ OC 



G(k,icj) 



as(itj) 
du! 



= 



which is necessary to provide the Luttinger theorem to be true. 

We use the following ansatz for the self energy is in the matrix (to, to ) form: 



E(icj) 



(iu) + A* £(2^ (iuj) * [l - B{iu) * S^^' (iuj) 



(13) 



where the second-order self energy S(^^(zcj) is defined in terms of the bath Green function Qo^iuj) in according with 
Eq.| with W = (see also Eq.|): 

In spirit of the approach p5| the ^-matrix should be defined to provide the exact high-energy (lu — s- oo) limit of 
T,{iu!). The best way to receive such asymptotic is the using of the equations of motion for the double-time retarded 
Green function with the analytical continuation on the Matsubara frequencies ||2^,Q. One has exact at ^ oo 
expression 
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where 

here [..., ...] and {..., ...} are the symbols for commutator and anticommutator, correspondingly, Hint is the Hubbard 
(interaction) part of the Hamiltonian. Note that in the multiband case the average Nmm'a contains the products of 
four electron operators and cannot be found exactly. Decoupling of these four-fermion averages according to the Wick 
theorem and comparing the result with the asymptotic of E^'^-': 



we obtain the following expression 



A = N * N° ^ 



where iVo-niatrix defined in the spin-polarized case as: 



{rrii} a' 

and in the non-magnetic case it simplifies to: 

The expression for A^-matrix differing from that for by the replacement of the occupation matrix 'nP n. Note 
that matrix A appears to be non-Hermitian. In the non-degenerate case this expression appears to be exact (see [^) 
due to the identity 



It can be quite accurate also in the general multiband case. 

Coefficient matrix B is designed to fix the exact atomic limit of the interaction self energy, Eq.|l^. There are 
other problems with coefficient B in the multiband case |^^. While in single band model one can find an analytical 
expression for the constant B | p5| ], in the multiband case this parameter should be cj-dependent, owing to the frequency 
dependence of the atomic self-energy, Eq.|l^. We decide to find numerically the non-Hermitian matrix B{iuj) from 
atomic limit of Eq.|l^ using the exact Sat(ia;) with a constraint for the density matrix n. In this limit Y,^'^\iuj) in the 
non-magnetic case has the form: 

y(2)at , . X _ jjd TT {frriji^ ~ fm2 ~ /mi) + /mi/m2] 

f ^ ^ f^O ^m2 ' ^ma ^mi 

here and e„i- are diagonal occupation numbers and energies of /i"* . In this case we have: 

B{iuj) = [E(2)°*(ic^)l - rE°*(iw) - gff^^j * A (14) 

As a simple example for such scheme we compare on Fig.|l| DMFT to exact-diagonalization for the Anderson model 
of two-site, two-bands with one correlated site U — 4, Sf — —4 and one "free-site" with Eq = and hybridization 
between the sites V = 0.25 [Q. For convenient we assume that all parameters for our model calculations are in "eV" 
energy units. The corresponding Hamiltonian for Anderson impurity model has the following form: 

Himp = £/ ^ ] fmafma + V ^ ^ {fmcr'^mcT + C^^fma) + U ^ ^ f fm] f ,-filf™-[ 
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It is not a problem to find an exact Green function for this model (non-symmetrize many-body Hamiltonian has the 
dimension 256x256) and compare with approximate calculations. We see that the agreement between exact solution 
and our DMFT results is quite good even for a large filling {ntot > 1; in this case ny = 0.76, ntot ~ 2 ). Also note that 
the atomic Green-function in Fig.|l| for correlated site has the "three-peak" structure for this occupations (there are 
general 8 poles in Green-functions for two-band case) and not the two-peak structure as in the one-band model. The 
use of the numerical atomic Green fmiction for the i3(iLj)-matrix calculation is quite important even for qualitative 
agreement with exact results for such model at the filling larger than one electron per site 



III. RESULTS FOR TWO-BAND MODEL. 



In this section we compare the three different LDA-I--I- approaches described above for a two-band system. We 
used the simplified two-dimensional model for High-Tc superconductors for d^2_y2 = x and d^^ = z orbitals 
If one can skip z-orbital it will be standard single-band nearest neighbor hoping ("t")-model. LDA band structure 
calculation for High-Tc materials shows the large contribution of Cu d^^ orbital to states near the Fermi level p9[ . 
Therefore the situation with two correlated valence bands could be possible in this materials. Although we knew 
that for the realistic description of Cu dx2_y2 state in single band model one need to include next-nearest hoppings 
(coming from interactions with 02p- and Cu4s-orbitals [[49[ ) we used here the simplified tight-binding model for two 
correlated bands "x" and "z" within nearest-neighbour hopping approximation. The one electron Hamiltonian has 
the following form [pij : 



h{k) = 



-2txx{cos kx + cos fcj,) + A, —2txz{cos — cos ky) 

-2t, ' - 



-2txz{cos kx — cos k. 



txzicos kx — COS ky) 
tzz{cos kx + COS ky) 



The hopping parameters are related via simple Slater-Koster ratio: txx — 1, tzz — 0.3, txz — 0.4. Again we assume 
that all TB-parameters are in "eV" energy units, while the value txx '~0.5 eV would be more realistic |Q. It is 
important to take into account the energy shifting parameter A since Cu dx2-y2 bands located higher than Cu dz2- 
one, so we use A = 4. For the Coulomb energy our parameterization corresponds to the following matrix elements 

{tTIx 7^ ^2)' ^mimimimi U -\- f^mim2mim2 ^ i Umim2m2mi (/,and t^mim2mimi 1^ this CaSe the 

symmetrized bare vertices has the following form (the basis function numbering as: xx, xz, zx, zz): 



U' = 



U + J 5 J 5 J 2U-J- 

5 J 2J -U SJ 
5J 2J~U 6 J 

2U~J 6 J 6 J U + J 

/U + J SJ SJ 

SJ UJ + u) UJ + u) 
SJ liJ + u) liJ + u) 

SJ 



V 



SJ 







/o 


SJ 







SJ 





u + j) 







-U-J ~SJ 

-SJ 

-SJ -U 

-J -SJ 






i 



(J-u) UJ-u) 











We investigate this model for different U parameters: U — 2 — 8 with the fixed values of J = 0.5 and SJ = 0.1. The 
total number of electrons are ntot = 1-4 which approximately corresponds to fully occupied z-bands and almost half- 
field a;-band with 10% of holes. We use the 32 x 32 mesh for the summation over the Brillouin zone and 4000 — 8000 
Matsubara frequencies with the cutoff energy equals to 20 — 40 times the bandwidth. On the Fig.|| we show the 
results of self-consistent FLEX two-band calculations for U = 2,4,6 and 8. Density of state (DOS) was obtained 
from the Green function, extrapolating Matsubara frequencies with Fade approximation to the real axis. Note 
that the bare bandwidth for txz = is equal to 8t in the case of two dimensional square lattice and corresponds 
to 8 and 2.4 eV for x and z bands. One can see the narrowing the DOS-peak near Fermi energy (Ep) for a;-band 
and boarding the total x- and z- subbands as U increased. An interesting feature of this two-band model is seen 
in moving the peak from occupied z-bandstowards the Fermi level with increasing correlation strength U, and a 
pinning of the van-Hove peak from the x-band just above Fermi level. This drastic change of the x- and z- band- 
shape maximize the particle-hole interband susceptibility (inversely proportional to the energy distance between band 
peaks) and therefore increase the fluctuation contribution to FLEX-energy. The value of the average inverse mass 
enhancement factor Z = [1 — dT, (e) /de]J^Q is equal to 0.45 and 0.83 for x and z orbitals for U=4. Spectral function 
A(k, e) — — l/7rTr/mG(k, e) for three different k-directions in two-dimensional Brillouin zone is shown in Fig.^. One 
could see the renormalized dispersion of two bands: x-band from approximately -5 eV at X-point to 5 eV at M-point 
and z-band from approximately -8 eV at X-point to 3 eV at M-point. Similar to all High-Tc models there is an 
extended van-Hove singularity in x-bands just at the Fermi energy near X-point. 
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The results of self-consistent DMFT calculations for U = 4 and 8 are presents in Fig.^. In this case we use only 
first seven Matsubara frequencies in the Eq.^ and the constant B-matrix for the rest frequencies. One can clearly 
see some differences to the corresponding FLEX results, which is related to a sharpness of the DOS near Ep and 
more pronounced three peak structure of partial x-band DOS for U=8. We plot also the DOS corresponding to 
the atomic Green function for U=8 with the four poles near to Fermi energy out of eight poles in the paramagnetic 
two-orbital atom. Corresponding spectral function for the same directions in Brillouin zone is presents in the Fig.^ 
for U=8. There is a sharp quasiparticle dispersion near the Fermi level and broad incoherent background above Ep 
at the energy of about 7 eV near M-point. Extended van-Hove singularity at X-point become more pronounced. We 
plot also the momentum distribution function n(k) in (1,1) direction in the Fig.^. From the quasiparticle dispersion 
of the x-band along T — AI direction (FigJ|), we expect the Fermi surface crossing almost exactly at the half-way 
between this two points. The momentum distribution function (Fig.^ just confirm this situations and shows that the 
Fermi step (our simulation temperature T=0.06t) is smaller than one and agree with the calculated value of the mass 
renormalization factor Z—0A3 for a;-band. 

Finally the HIA-solution for this two-band model for [/ = 8 is shown in Figj^. In this case we have the dielectric 
DOS with narrow atomic-like resonances. It is interesting to mention that the structure of the atomic Green function 
in the DMFT approximation (Fig.^ quite close to the HIA-solution shifted approximately by 2 eV down. We would 
like to mention that the application of HIA-scheme is reasonable only for U>>W and not for U=W as in this case. 



IV. LDA++ CALCULATIONS FOR REAL SYSTEMS 



The self-consistent LDA-I--I- calculation for real systems pose a serious computational problem. One need to operate 
with the susceptibility matrix, which is of dimension x and depends on the Matsubara frequencies. For a 
illustrative purpose we have calculated the electronic structure of classical Mott-Hubbard insulator NiO and mixed 
valence 4f-compound TmSe in HIA scheme for LDA-I-+. We use non-selfconsistent HIA approximation with the 
simplest constrain for only diagonal £/ in Eq.^ to have rid for NiO or n/ electrons for TmSe from self-consistent 
paramagnetic LDA-calculations. The Coulomb matrix is expressed via effective Slater integrals: 

k 

where < k < 21 and 

An ^ 

ak{mi,m[, 1712,1712) = 2k + I ^ ^^'"^ ' "^'"^ ' lm[){lm2 \ I ^"4) 

q— — k 

We used the following effective Slater parameters, which define screened Coulomb interaction in d-shell for NiO: 
F° = 8.0 eV, = 8.2 eV^^ = 5.2 eV, and in f-shcU for TmSe: F° = 5.7 eV, = 9.1 eV, F^ = 5.7 eV, 
F^ — 4.7 eV, (see e.g. [ |l6| , pl|j5^ ) . We start from the non-magnetic LDA calculations in the LMTO nearly-orthogonal 
representation psf for experimental crystal structures of NiO and TmSe. The minimal basis set of s,p, c?-orbitals 
for NiO and s,p, d, f- orbitals for TmSe corresponds to 18x18 and 32x32 matrix of the LDA Hamiltonian /i(k). The 
occupation number for correlated electrons are 8.4 electrons in the d-shell of Ni and 12.6 electrons if the 4f-shell of 
Tm. Using the corresponding atomic self energy for Ni-atom and Tm-atom the total DOS for NiO and TmSe have 
been calculated from Eq.^. In Fig.^ we compare the paramagnetic LDA results with HIA LDA++ scheme. It is 
well known that paramagnetic LDA calculations can not produce the insulating gap in nickel oxide: the Fermi level 
located in the middle of the half-filled eg bands |jl^. In the HIA approximation to LDA-I--I- approach there is a gap 
(or pseudogap in the Fi^ due to temperature broadening) of the order of 3.5 eV even in this "nonmagnetic" state. 
This gap and the satellites at -5 and -8 eV are related to the structure of atomic Green function shown in the lower 
panel of the Fig.^. 

In Fig.Hwe compare the calculated DOS for TmSc with experimental XPS spectrum ||5^. The HIA approximation in 
this case well reproduces the ladder-type photoemission spectrum, comes mainly from combinations of two multiplets 
structure of Tm^+(/^'^) and Tm'^+(/^^). This example demonstrate how the LDA-I--I- scheme can combine the many- 
body atomic physics with band structure methods. Normal LDA-band structure for rare-earth systems corresponds 
to the narrow f-peak at the Fermi level and could not describe the experimental XPS-spectrum which has the "f"- 
resonances over wide energy range of the order of 12 eV. At the same time, HIA is not adequate to describe correctly 
the fine features of the electron structure near the Fermi level. It is known that TmSe is really the narrow- 
gap semiconductor. According to the most developed model approach to mixed-valence semiconductors [pSP] the 
appearance of this energy gap is caused by both hybridization and exciton effects due to the Coulomb attraction of 5d 
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conduction electron and 4f hole. This effect cannot be described in HIA approximation. Nevertheless we believe that 
the description of the electronic structure of f compounds including mixed-valence ones on the large energy scales is 
important itself and in this sense the results presented here demonstrate the usefulness of LDA-I-+ approach for the 
description of real strongly correlated systems. 



V. SUMMARY 



We have formulate a general LDA++ scheme which takes into account dynamical electron fluctuations in the case 
of correlated d- or f-states. The most accurate approach is the DMFT-band structure method, while more simple 
FLEX and HIA scheme can be as well useful for investigation of correlation effects in real systems, in the cases of 
rather weak and rather strong interaction respectively. In principle one could combine the idea of the bath-Green 
function in DMFT-scheme with the simple expression for the self-energy in the FLEX approximation. In this case 
S = 'E^^^^[Qo] and we expected effectively reduce the effects of vertex-corrections in the FLEX scheme. 

Here we compare the LDA+-I- approach with more simple LDA-I-U one. First of all, to describe Mott insulators 
in LDA-I-U approach (as well as in SIC approach) it is necessary to assume magnetic and (or) orbital ordering [p^ . 
In LDA-I--I- it is possible to consider the paramagnetic Mott insulators in the framework of ab initio calculations. 
Moreover, it is possible to obtain not only the Mott- Hubbard gap in the electron spectrum but also satellites and 
multiplet structure (see e.g. the results for TmSe and NiO in the previous Section). 

The correlation effects results from the frequency dependence of the self energy ( the non-Fermi-step form of the 
distribution function for quasiparticles, the mass enhancement, the appearance of many-electron Kondo resonanses etc) 
can be obtained and investigated in LDA-|--t- approach but not in LDA-I-U. Our results for two-band model provide the 
interesting examples of such behavior. In particular it is worthwhile to note such features as the narrowing of the van 
Hove singularity and its "pinning" to the Fermi level (which is important for the physics of High Tc superconductors 
p6| and can be described already in multiband FLEX approximation), the "three-peak" structure of the spectrum in 
the vicinity of Mott insulators ("Kondo resonance" and mid-gap states which can be described in DMFT approach). 

We hope that the approximations described here may be useful for the ab initio calculations of electron structure 
of a great variety of strongly correlated electron systems including doped Mott insulators, rare-earth metals and their 
compounds (in particular mixed- valence ones), high-temperature superconductors and many others. 
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VII. APPENDIX A 



Here we present the proof of the equivalence of a " cavity method" and impurity problem in DMFT for multiband 
case. We start with the expression for the Green function matrix on the zero site in the cavity method (see 
Eq.(35)) 

[Go{i(^)]~^ ^ iu! + - hat - R, 

where 



with hat being the one-electron part of the intra-atomic Hamiltonian. All this functions are the matrix in (m, m') 
indices as well as diagonal matrix in spin-space. Note that G,-^'' is the Green function between sites i and j on the 
lattice with the site zero being eliminated 

G\j^ = Gij — Gi^GQ^Goj (16) 
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Using the Fourier expansion of all the quantities over the Brillouin zone and substituting Eq.fiq to Eq.na one has 



where 



and 



At the same time (see 0) 



G(k) = [A-t(k)]-i 



A ^ iuj + fi — hat — S(ici;) (17) 
Taking into account that J2^0^) = '^^'^ obtains 

k 

L = Y, Wk) - A + A] [A - i(k)]"' = -1 + AGqo 

k 

and 

= -l + GooA 

One can obtain by the similar way the result for M-matrix: M = AL^ 
Substituting all this formulas to Eq.|l5| we have 

R^iA + Goo' - A)L^ = -Goo' + ^ 

using Eq.|l^, we can finally write: 

[^o(iw)]~' = iuj + fj, - hat - A + Goo' = S(icj) + Gqo' 

It means that C/o(*^) is the Green function of the impurity problem with the on-site one-electron Hamiltonian 
hat + S(ict;) on each "non-zero" site and hat for zero site. 
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FIG. 1. Energy spectrum for two-band two-site Anderson model in exact diagonalization and DMFT scheme as well as 
atomic green function for correlated cite. 

FIG. 2. Density of states for two band model in the FLEX scheme for different U-values. Full and dashed lines indicate 
partial DOS for x- and z- orbitals. 

FIG. 3. Spectral function (the two band FLEX model, U=4) for three different directions in the two-dimensional square 
Brillouin zone: F = (0,0), X = (0.5,0), and M = (0.5,0.5) in units of 27r/a. 

FIG. 4. Density of states for two band model in the DMFT scheme for different U-values as well as atomic Green function 
for U=8. Full and dashed lines indicate partial DOS for x- and z- orbitals. 
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FIG. 5. Spectral function (the two band DMFT model U=8) for three different direction in the two-dimensional square 
Brillouin zone 



FIG. 6. The momentum distribution function n(k) in (1,1) direction for the two band DMFT model U=8. 

FIG. 7. Density of states for two band model in the HIA scheme for U=8. Full and dashed lines indicate partial DOS for x 
and z orbitals. 

FIG. 8. Density of states for paramagnetic nickel oxide in the LDA and HIA approximations as well as Ni-atom Green 
function. 

FIG. 9. Density of states for TmSe in HIA sheme in comparisson with experimenral XPS-spectrum and results of 
paramagnetic LDA-calculations. 
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